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Abstract 

The cross section for elastic vector meson production in photon- 
proton interactions at large t is considered using the exact analytic 
solution of the BFKL equation in the azimuthally symmetric n = 
limit. We use a non-relativistic model for the vector meson production 
and find a small shrinkage in the i-distribution with increasing energy. 
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1 Introduction 



We study the process 7P — > V + X (where V is a vector meson and X 
denotes the products of the proton dissociation) within the Regge hmit 
assuming a large moment transfer, s ^ —t ^ Qq, i.e. large angle meson 
production at high centre-of-mass energies. Qo denotes the typical hadronic 
scale of ~ 1 GeV. The largeness of —t allows us to apply the solutions 
to the nonforward BFKL-equation which are given as a set of conformally 
invariant eigenfunctions |l|] (BFKL-Pomeron). This process has been studied 
previously in ||2| , where an analytic solution of the BFKL equation was given 
for asymptotically large s and for —t either very large or very small (relative 
to the larger of the meson mass and photon virtuality). Comparisons with 
an iterative numerical solution were also made. An appreciable cross section 
was found and the prospect of revealing the BFKL dynamics at HERA was 
established. 

Here we derive a formula for the cross section which is based on the 
solution of the nonforward BFKL equation in the azimuthally symmetric 
approximation. We present a detailed numerical evaluation of this formula. 

Throughout we work in the approximation that the light-cone wave- 
function of the vector meson is a (5-function, which partitions the light-cone 
momenta and the transverse momenta of the quark-antiquark pair (from the 
photon) equally. This should be appropriate for heavy vector mesons (e.g. 
J/^' and T) and possibly also for the lighter mesons (p, uj, </>), however, 
the choice of the wave function is not uniquely determined. One could also 
think of a softer distribution than a (5-function. 

For the coupling of the Pomeron to the proton we use the Mueller- Tang 
prescription where the pomeron couples to a single parton only. This 
is justified Q since, in the case considered here, the size of the Pomeron 
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(l/-v/^) is smaller than the size of the proton (l/Qo)- If the size of the 
pomeron becomes larger the full triple pomeron vertex |Q, ^ has to be taken 
into account. 

The main concern in this paper is the exact convolution of the Lipatov 
eigenfunction with the meson wave function. At large enough energies the 
azimuthally symmetric solution (n = 0) gives the dominant contribution 
and we can neglect eigenfunctions with n 7^ 0. We study the resulting t- 
dependence and look for a possible shrinkage in the t-shape with increasing 
energy. Comparison will be made with the previous analytic calculations 
and numerical studies |0. 



2 The exact BFKL-solution 

Let us begin by defining the cross section for the diffractive scattering of a 
virtual photon off a proton: 

dcrjjq Vq) ^^vr 2 
It = ^T^'-^^"'")' • 

Since the coupling of the BFKL-Pomeron to gluons is the same as the cou- 
pling to quarks except for a constant colour factor (81/16), we have intro- 
duced, as our basic subprocess, the scattering off quarks only. We have 
chosen to factorise the l/f^ behaviour since it allows us to write the remain- 
der of the cross section in terms of the scaling variable 

r=-t/{Q^ + ml.). 

is the photon virtuality (it may be zero) and my is the vector meson mass. 
The rapidity gap between the produced vector meson and the scattered 
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quark is proportional to the 'energy' variable z: 

z = In 



27r -t + Q2 + ' 
where s is the Mandelstam variable representing the centre-of-mass en- 
ergy (s = x'W'^ with W being the total hadronic energy). From the work 
of Lipatov |T| , we have 



du- 



where 
and 



(27r)3 {u^ + 1/4)2 

X{u) = -4(7s + 7^e V(l/2 + iu)) 



(2) 
(3) 



V{k,q) I (fpicfp2 



PipI 



l/2+iu 



Jk-pi+i{q-k)-p2 ^ ^^-^ 



(2vr)2 

The momentum transfer —t = q'^, and for the pomeron- meson vertex we use 



the non-relativistic form factor 

C ( 



V{k,q) 



1 



1 



2 \qf + {k-q/2f qi^q^l^ 



(5) 



where 



'Scleral 



and 4(7^ = + m-v- 



In this approximation, there is no helicity-flip amplitude and the rate of 
production of (longitudinally polarised) mesons off longitudinal photons is 
simply Q"^ /my times the transverse cross section. 

The factor {q) is computed using the subtraction prescription of Mueller 
and Tang i.e. 



(pk 

(27r)2 



Cppi(fp2 



\ l/2+iu / 1 \ 



-)■ 



Jk-pi+i{g-k)-p2 



(6) 
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It is straightforward to evaluate eq.(6): 

We now turn to the more comphcated evaluation of eq.(4). The A;- 
integral can be performed exactly (it gives a McDonald function). Subse- 
quently, it is convenient to use the variables R= {pi + p2)/'2. and p = pi — P2- 
After a Feynman paramctrisation: 



C ,2j, (\ p^+^-'-i^o(g||p)[x(l-x)]-V^+^- 

[x{R - p/2f + (1 - x){R + p/2)2]i+2^- 



By changing variables to i?' = i? — (a; — l/2)p we are able to perform the R' 
integral and the p angular integration. We get 

T^(l/2 + iu) \2j Jo 

roo 

X dp p2 Jo{qp{l/2 - x)) Ko{q\\p) K2Uqp,/[x{l - x)]). (9) 



The trick now is to use the following representation of the 5-function: 

"1/2+400 

/1/2-ioo 



/•1/2+ioo J„ 

s{p'-p")= f.p-'-''p"'-\ (10) 

Jl/2-ioo 2m 



to give 



"'^^^^ = ^ r2(i/T+i.) (1) / vKi-x)] /'^/l^ 

X pi-2Vo(gp(l/2 - x)) X2i.(W[a;(l - a;)]) 

X j dp'p'^'Ko{p'q\\). (11) 

After performing the p' integral, we are able to isolate the important q/q\\ 
factor. Putting ^ = qp we have 

ds f dx 2^' ( q\ 7rr2(l/2 + s) ( q^ 



9'^ /gy+'^ 7rr2(l/2 + .) U^' 



j e'-'Vo(C(l/2 - x)) - x)]). (12) 
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The ^ integral can be done to give 



tV ( \ p ^ f f 



TT f as f ax 1 q \ Q \ 



r2(l/2 + ii/)7 2m J [x{l - x)f/^^^ \q\\ J \4. J 

X r2(i/2 + s)r(i-s-iz/)r(i-s + ?;z^) 

X F{l-s-iu,l-s + iu;l;-^-^). (13) 

4x(l — X) 

Changing variables to y = (1 — 2x)'^ allows the x integral to be performed, 
giving 

rV.. ^ VT r(l/2)r(l/2 - z^) ( q^Y" f ds ( qV^'' 

''q3r^{l/2 + iu)r{l-iu)\4 J J 2m \qj 

X r^{l/2 + s)r{l-s-iu)r{l-s + iu) 

X ^F2{s-ii',l- s-ii',l/2;l-ii';l). (14) 
Using Watson's theorem, this can be simplified to the important result: 



(15) 



q^r{l/2 + ii^) \ A J ii/2-ioo 2m 
r(l - s - iv) 
^ r(s/2 + ii^/2 + l/2)r(s/2 - iiy/2 + 1/2) 

r{l- s + iu) 
^ r(l - s/2 + iu/2)T{l -s/2- iv/2) ' 

Putting equations (7) and (15) into (2) we arrive at our final answer, 

/CO ,,2 /•1/2+joo fla 

-oo (l/2 + l/4)^ A/2-ioo 27rz 

r2(i/2 + s)r(i -s-iz^) 



r(s/2 + iv/2 + l/2)r(s/2 - iz//2 + 1/2) 
T{l-s + iu) 



" r(l-s/2 + zi//2)r(l-s/2-zi//2)' ^^^"^ 
In fig.l, the solid lines show the result of performing the s and v integrals 
exactly, i.e. they correspond to the exact solution to the BFKL equation for 
n = 0. The r dependence at different z values is shown. 
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3 The limits r > 1 and r <C 1 



The s-plane integration in eq.(16) is quite straightforward to compute. It 
leads to a power series in r (or 1/r). For r > 1, we must close the contour 
in the left half plane and for r 3> 1 it is the double pole at s = — 1/2 which 
provides the leading power behaviour. Similarly, for r < 1, we must close 

in the right half plane and the leading behaviour is determined by the two 
poles at s = 1 lb izv. 
For r 3> 1, we have 



The V integrals can easily be performed numerically. 

Fig.l illustrates how remarkably well these approximate forms reproduce 
the full solution (solid lines). For r < 0.3, cq.(18) reproduces the full solution 
over the whole z range (dashed line). For large r the agreement is excellent 
for large enough z - even for r as low as 0.5 for z = 0.8 (dash-dot line)! 
The vanishing of the cross section at t ^ 10.1 for z = 0.2 reflects a large 
negative contribution to the amplitude from large u (where the integral is 
very slowly convergent) and indicates the break down of the approximation. 

Prom eq.(17) we conclude that, at large enough r, the dependence upon 
r and z factorizes, i.e. the t-dependence is energy independent in this case. 
The situation, however, changes for r < 1. Eq.(18) is dominated by the 




(17) 



and for r <C 1 




(18) 
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exponent: 

ex.p{x{i^)z + (3/2±zi/) Inr} (19) 

We can make a simple estimate using the saddle point method with two 
different approximations of the ^-function. In order that this exponent be 
dominant, we require Inl/r and z to be large enough. 

In the first case we assume In l/r to be not too large compared to z and 
expand x{^) arround u = 0: x(i^) ^ 8 In 2 — 28^(3)i/^. The saddle point is 
then 

56C(3) z 

This limit refers to the BFKL solution at large z. The formal limit of 
z — DO will be discussed in the next section. 'Large enough z' means that 
z is much larger than ln(l/r)/[56C(3)] so that the saddle point lies in the 
region (i/ 0). The saddle point integration then leads to 

82ln2 / \ 

^(,,,)._exp^-^j^]. (21) 

The second case wc consider is for In l/r S> 2;. In this limit, we need to 
expand the x function around ly ~ zbz/2: x(^) ~ 2/(1/2 zb iv). The saddle 
point now reads: 



The solution for the amplitude is then the well known double leading log 
result, i.e. 

^2 exp{2V2^1n(l/r)} 
>(2;,T)~r ^ . {26} 

This result has an intuitive interpretation. For small —t and large + my 
the dominant contribution comes from strongly ordered transverse momenta 
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within the BFKL-Pomeron, starting at —t and ending at + mY. As soon 
as the internal transverse momenta are larger than —t the effect of t in 
the loop integrations is negligible and the solution is identical to the t = 
evolution with —t entering only as the starting scale. The cross section scales 
as expected, i.e. because of the factor it behaves like C^/((5^ + rriy)'^. 
Thus it is strongly suppressed at large (it is a higher twist contribution) . 

In order to estimate the magnitude of the shrinkage we calculate an 
effective slope Oi'^jf following the soft Pomeron approach by differentiating 
the logarithm of the amplitude with respect to t and Ins, i.e. 

/ 1 d"^ , ^, . , 



, _ 1 oas mi/r 



«e// = ^ — . /_^ • (26) 



When In 1/r ^ 56C(3)z we can use Eq.(21). In which case, 

1 3as In 1/r 
|i|^56z2C(3)' 

When In 1/r 3> 8z we can use Eq.(23). In which case, 

1 SOs 1 
|t| 27r ^2 zln(l/r)' 

A different approach , based on the diffusion in the impact parameter 
space leads to a dependence of a'^jj on and Ins similar to Eq.(26). (i.e. 
a'^jf ~ y/os/lns). 

In fig. 2 we plot the energy dependence of the dimensionless quantity, 
— d'^ In J^/dzdr. The effective slope is obtained simply by multiplying by 
3Q!sr/(27r|t|). The solid line corresponds to the slope extracted from Eq.(18). 
The dotted line corresponds to the effective slope given by Eq.(25) and the 
dashed line is given by Eq.(26). From the figure, we can see clearly that the 
effective slope is small for all accessable values of r and z, i.e. Og// ~ 
GeV~^ for —t>l GeV^ and z > 0.1. Recall that the slope of the soft 
pomeron trajectory advocated by Donnachie & Landshoff is ~ 1/4 GeV~^ 
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[^. Also we note that the saddle point evaluations do not lead to very 



good agreement with the full solution for sensible values of z and r (we 
have checked that the saddle point evaluations do eventually agree with 
the numerical evaluation for large enough z (eq.(25)) or small enough r 



To illustrate the shrinkage in a more familiar way, in fig. 3 we show the 
actual i-distribution for J/*!' photoproduction (Q^ = 0). The distribution is 
shown for the 79 Vq process at different jq energies, s. For comparison 
the cross sections at each energy are normalised to their value at —t = 1 
GeV2. We took = 1/4 and s = where W = 50 GeV (solid line), 100 

GeV (broken line) and 200 GeV (dotted line). These energies are typical of 
the HERA kinematics. The corresponding 75 cross sections at —t = 1 GeV^ 
are: 16 nb/GeV^, 72 nb/GeV^ and 360 nb/GeV^ respectively. The 7p cross 
sections are (in the valence quark region of large x' ~ 1/3) approximated 
by the individual quark cross sections multiplied by a factor of 3. Note that 
the z range probed by these energies is 0.5 ^ z < 0.9. 

4 The limit 2; > 1 at fixed r 

The asymptotic BFKL limit is that of z ^ 1 at fixed r. In the region where 
eq.(20) is valid we can present fully analytic formulae for large and small r, 
by recognising that the dominant contribution to the z^-integral only arises 
from the region around = 0. As already detailed above, in this limit, 
X{u) 81n2 - 28C(3)z^2 ^nd so 



(eq.(26))). 




(27) 



Thus we obtain, for r » 1: 



16^3/2 

r4(i/4) 



In 16r 



[7zC(3)]3/2 



(28) 
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This is as in Eq.(21) after assuming z to be large enough to neglect any 
z-dependence of the t behaviour. 
For r < 1: 

gSz In 2 

T{z, t)^C^ [In 1/r + 61n2 - 4]r3/2^_-^-^. (29) 

Eqs. (28,29) were also found in ref.[Q]. The z ^ 1 limit is represented by 
the dotted line in fig.l, and can be seen to provide a faithful representation 
of the exact solution for z ^ 0.8. We should point out that the dotted lines in 
the figure actually correspond to the numerical evaluation of the amplitude 
in the v = ^ limit, i.e. they are valid for all r. However, outside of the 
region 0.2 < r < 1 eqs. (28,29) produce essentially the same results. 



5 Conclusions 

We have shown that the exact analysis of the BFKL-pomeron in the az- 
imuthally symmetric approximation (n=0) leads to a small shrinkage of the 
t-dependence at small r ( < 0.1) which dies away with increasing z. For 
z ^ 1 the asymptotic formulae of eqs. (28,29) start to become reasonable. 

At small z < 0.5 and small r < 0.5, the exact n = solution is in nice 
agreement with the results of a fixed order BFKL summation, represented 
by the diamonds in fig.l. The diamonds were first presented in ref.Q and 
correspond to computing numerically the BFKL corrections to the Born 
(2-gluon exchange) amplitude up to and including terms ~ z^ . For larger 
z, it was shown in ref.Q that the perturbation series is slowly convergent. 
This is why the diamonds lie below the full n = solution. Certainly 
the numerical results of ref.j^ strongly suggest that, for r < 0.5, the n = 
solution to the BFKL equation derived here is an excellent approximation to 
the full solution, even in the region of small z (where the Born amplitude is 
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sufficient). The region of small r is experimentally accessable in light meson 
electroproduction and in heavy meson photo- and electro-production. 

However, the situation changes when r is increased. The diamonds in 
fig.l no longer agree with the n = solution (solid line). In ref.[Q] it was 
found that the large r region is difficult to compute, and some progress was 
made in generating reliable results. One of the main difficulties arises due to 
the presence of a dip in the r-distribution [Q]. The origin of this dip can be 
traced back to the peculiar behaviour of the Born amplitude, i.e. it vanishes 
for r = 1. However, no dip is present in the n = solution. When z is 
increased the dip moves towards larger r (but does not seem to disappear). 
More work is needed to understand the behaviour of the amplitude in the 
region of the dip. Experimentally, the region of the dip is accessable in the 
high-t photoproduction of light vector mesons. 

It is of great interest to extend the present calculation towards smaller —t 
(^ Qo)- This kind of extension requires one to go beyond the Mueller- Tang 
approach, pursuing the improvements made in Q. We expect the shrinkage 
in the t-dependence to become stronger ||5|. 

The first HERA data on the the diffractive photoproduction of J/^ 
mesons contains a number of events at high (i.e. ^ 1 GeV) The 
data sample is too small at present to draw any strong conclusions and we 
look forward to the increase in statistics. 
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Figure Captions 



Figure 1: The solid line represents the full BFKL-solution (eq.(16)), the 
dotted line shows the numerical evaluation of eq.(16) about the ~ saddle 
point. The dashed line corresponds to eq.(18) (i.e. r ^ 1) and the dash- 
dotted line to eq.(17) (r S> 1). The diamonds are from ref.Q and correspond 
to a numerical calculation of BFKL corrections to the Born contribution. 

Figure 2: The dimensionless ct'^fj from eq.(24). The solid line shows the 
numerical evalution based on eq.(18), the dotted and dashed lines the cor- 
responding saddle point solutions of eqs.(25) and (26). 

Figure 3: The cross section normalized by its value at \t\ = 1 GeV^ for 
three different energies: W= 50 GeV (solid line), 100 GeV (broken line) and 
200 GeV (dotted line). 
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